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Batteries  are  being  used  in  increasingly  complicated  configurations  with  very  demanding  duty  schedules. 
Such  usage  makes  the  use  of  batteries  in  multi-cell  configurations  to  meet  voltage,  power,  and  energy 
demands  in  a  very  stressful  manner.  Thus,  effective  management  and  control  of  a  battery  system  to  allow 
efficient,  reliable,  and  safe  operation  becomes  vital,  and  diagnostic  and  prognostic  tools  are  essential.  Yet, 
developing  these  tools  in  practical  applications  is  new  to  the  industry,  difficult  and  challenging.  Here  we 
present  a  novel  mechanistic  model  that  can  enable  battery  diagnosis  and  prognosis.  The  model  can 
simulate  various  “what-if’  scenarios  of  battery  degradation  modes  via  a  synthetic  approach  based  on 
specific  electrode  behavior  with  proper  adjustment  of  the  loading  ratio  and  the  extent  of  degradation  in 
and  between  the  two  electrodes.  This  approach  is  very  different  from  the  conventional  empirical  ones 
that  correlate  the  cell  parameters  (such  as  impedance  increases)  with  degradation  in  capacity  or  power 
fade  to  predict  performance  and  life.  This  approach,  with  mechanistic  understanding  of  battery  degra¬ 
dation  processes  and  failure  mechanisms,  offers  unique  high-fidelity  simulation  to  address  path 
dependence  of  the  battery  degradation. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Battery  performance  and  its  degradation  have  significant 
consequence  on  reliability  and  safety  of  a  battery  system’s  function, 
which  shall  determine  its  feasibility  and  acceptance  in  the  electric 
drive  vehicle  and  electrical  energy  storage  applications.  To  under¬ 
stand  battery  degradation  mechanisms  in  a  clear  cause-and- 
consequence  context  remains  very  challenging  today,  and  such 
a  long  struggle  in  the  research  community  is  likely  due  to  the  lack  of 
simple  and  effective  diagnostic  tools  to  characterize  degradation  in 
situ  through  aging.  Many  postmortem  analyses  might  reveal  the 
cause  of  degradation  and  failure  of  a  battery,  but  they  are  not  useful 
in  providing  temporal  resolution  for  practical  battery  monitoring, 
protection,  and  diagnosis  during  operation,  nor  can  they  provide 
useful  information  for  prognosis.  In  our  view,  an  effective  battery 
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diagnosis  needs  two  major  elements:  a  reliable  method  to  deter¬ 
mine  the  state  of  the  battery  with  sufficient  temporal  resolution 
and  a  high-fidelity  model  to  predict  the  performance  under  a  duty 
regime.  To  achieve  such  a  goal,  here  we  present  a  unique  mecha¬ 
nistic  model  that  can  reveal  specific  battery  chemistry  and  its 
degradation  modes.  A  degradation  mode  includes  a  degradation 
mechanism  and  consequences  on  battery  performance,  such  as 
degrees  in  the  fade  of  capacity,  rate  capability,  and  peak  power;  and 
cell  properties,  such  as  the  increase  in  polarization  resistance;  with 
sufficient  temporal  resolution  to  tell  the  differences  under  various 
operating  conditions. 

The  most  effective  and  simple  method  of  monitoring  battery 
behavior  is  through  the  measurements  of  battery  voltage,  current, 
and  temperature.  In  recent  years,  we,  among  others,  noticed  that 
the  voltage  and  capacity  responses  of  a  Li-ion  cell  in  a  duty  regime 
could  be  analyzed  to  derive  information  that  can  be  used  to  esti¬ 
mate  the  state  of  a  battery  and  its  degradation  [1-12].  Yet,  most 
studies  still  focus  only  on  capacity  loss  against  cycle  aging.  Cell 
voltage  not  only  depicts  what  electrochemical  reactions  are 
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occurring  in  the  cell,  its  variations  also  reflect  how  the  reactions 
progress  in  a  duty  regime.  The  ability  to  analyzing  such  cell  voltage 
responses  in  different  duty  regimes  is  the  key  to  understand  cell 
behavior  and  degradation.  To  do  so,  a  unique  capability  that  we 
advocated  recently  is  the  analysis  of  incremental  capacity  (IC) 
[1-6],  i.e.  dQJdV  =  /(V),  in  diagnosis  on  a  variety  of  commercial  Li- 
ion  batteries,  including  those  made  of  LixNio.8Coo.15Alo.05O2  (NCA), 
LixFeP04  (LFP),  LixNii/3Mni/3Coi/302  and  LixMn204  composite 
(NMC  +  LMO),  as  reported  in  our  previous  work  [1-6].  Others,  like 
Honkura  et  al.  [7,8]  and  Bloom  et  al.  [9-12],  prefer  the  analysis  of 
differential  voltage  (DV),  i.e.  dV/dQ. =/(()),  and  find  it  useful  as  well. 
Although  both  methods  should  yield  similar  information  (in 
a  reciprocal  fashion),  there  is  a  distinct  difference  between  the  two: 
the  IC  refers  to  the  cell  voltage,  which  could  be  a  direct  and  primary 
indicator  of  the  battery  state;  whereas  the  DV  refers  to  the  cell 
capacity,  which  is  an  accumulative  and  secondary  indicator  that 
could  vary  with  aging  and  degradation  and  lose  its  reliability  as 
a  reference  in  the  course  of  aging. 

Disregarding  the  differences  in  the  two  techniques,  we  also  note 
that  the  interpretations  of  dQJdV  or  dV/dQ.  obtained  from  the  cell 
characterizations  are  not  trivial  at  all  to  reveal  cell  behavior  and  its 
degradation  modes.  Particularly,  since  both  electrodes  in  the  cell 
contribute  to  the  features  (peaks  and  their  shape)  in  the  IC  and  DV 
curves,  it  is  difficult  to  separate  the  electrode  behavior  and  its 
contribution  to  the  cell  degradation.  A  technique  that  can  resolve 
this  issue  would  be  highly  desirable  for  battery  monitoring,  control 
and  management.  In  our  earlier  work  we  solved  this  issue  by  using 
an  analytic  approach  with  qualitative  results  [5].  In  this  work,  we 
developed  a  synthetic,  mechanistic  model  using  half-cell  data  ob¬ 
tained  from  experiments  or  literature  to  simulate  the  evolution  of 
the  IC  and  DV  signatures,  so  we  can  quantify  electrode  contribu¬ 
tions  in  cell  degradation  modes.  This  approach  is  applicable  to  most 
of  the  known  chemistries  used  in  the  commercial  cells  and  it  can 
generate  various  “what  if’  scenarios  for  cell  degradation  modes  as 
case  studies  to  help  us  understand  how  cell  degrades  in  a  quanti¬ 
tative  manner. 

This  approach  conceptually  follows  the  one  described  by 
Christensen  and  Newman  in  the  early  2000s  [13,14]  and  used  by 
others  such  as  Zhang  et  al.  [15,16]  and  Delacourt  et  al.  [17,18]  to 
simulate  cell  degradation.  A  common  theme  in  these  approaches  is 
the  use  of  electrochemical  models  to  simulate  half-cell  electrode 
behavior  to  compose  the  full  cell  charge  and  discharge  curves.  DV 
or  IC  curves  could  be  derived  in  principle,  yet  results  of  such  deri¬ 
vations  have  not  been  clearly  reported  in  the  literature  to  date.  In 
contrast,  some  groups  used  experimental  half-cell  data  and  fitting 
techniques  to  compose  dV/dQ.  curves  [7,8,11,19-21]  for  some 
specific  studies.  These  experimental  approaches  may  provide  the 


benefits  for  case  study  and  analysis;  however,  the  applicability  of 
the  knowledge  to  other  chemistries,  cell  designs,  or  operating 
conditions  is  limited,  since  they  are  after  all  empirical  ones.  Often, 
when  substantial  changes  in  cell  chemistry,  design,  and  operating 
conditions  occur,  the  entire  process  has  to  be  re-invested  again  to 
receive  required  fidelity  and  accuracy,  which  is  quite  costly  and 
labor  intensive. 

The  approach  reported  in  this  work  combines  both  modeling 
and  experimental  techniques  to  provide  a  universal  tool  for  diag¬ 
nosis  and  prognosis.  Our  approach  employs  a  modified  equivalent 
circuit  model  (ECM)  from  which  the  cell  performance  is  emulated 
from  two  separated  half-cell  modules  each  constructed  from 
laboratory  experimental  data  [22-26].  This  approach  offers  the 
benefits  of  wide  applicability  of  the  model  to  various  cell  chemis¬ 
tries,  designs,  and  operating  modes;  as  well  as  the  high  fidelity  as 
inhered  from  detailed  extraction  of  the  experimental  data.  In  this 
work  we  use  the  LixFeP04  (LFP)  chemistry  as  a  case  study;  but  the 
approach  is  also  validated  with  data  obtained  from  our  prior  work 
on  other  chemistries  [1—6]. 


2.  Principles  and  model  descriptions 

The  mechanistic  model  proposed  in  this  work  is  composed  of 
two  layers:  a  top  layer  with  cell  configuration  and  inputs  from  cell 
operating  and  degradation  modes  and  a  sub-layer  with  half-cell 
modules  that  describe  the  electrode  behavior.  The  modules  use 
an  ECM  approach  similar  to  that  described  in  our  previous  work 
[22,23]  to  handle  the  electrode  chemistry  based  on  voltage  and 
capacity  variation  as  a  function  of  rate.  The  top  layer  uses  equations 
to  assemble  the  two  electrodes  into  a  full  cell  configuration  by 
considering  electrode  composition  and  loading  ratio  matching  as 
a  function  of  degree  of  degradation,  degradation  modes,  and 
operating  conditions. 

The  schematic  of  the  mechanistic  model  and  modeling  approach 
is  presented  in  Fig.  1.  In  contrast  to  conventional  electrochemical  or 
electrical  models,  in  which  the  degradation  results  are  outputs  of 
the  simulation,  ours  is  a  “backward  looking”  approach— the  “what 
if’  scenarios  of  degradation  are  the  inputs  to  the  simulation  of  the 
cell  behavior  in  aging.  The  architecture  of  the  model  accommodates 
electrode-specific  phenomena.  It  can  also  handle  complicated  cell 
and  electrode  configurations,  including  composite  electrodes  that 
comprise  multi-composition  mixtures  of  grains  or  core— shell 
micro-structured  grains.  The  model  is  written  in  MATLAB®  with 
architecture  configured  for  effective  functioning  in  computation, 
layer  integration,  and  flexibility  in  conducting  simulations  of  “what 
if’  scenarios. 
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Fig.  1.  Schematic  of  the  modeling  approach  in  this  work. 
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2.1  Half-cell  sub-layer  modules 


The  half-cell  module  in  the  model  is  a  key  element  in  this  work. 
For  a  typical  single-composition-based  electrode,  the  simulation  is 
performed  by  feeding  the  ECM  (Fig.  1)  with  parameters  derived 
from  half-cell  data  for  each  electrode  [22,23].  For  more  complicated 
electrode  architecture,  the  electrode  ECM  can  be  composed  from 
several  sub-ECMs  each  representing  an  active  material  composition 
in  series  (core— shell)  or  in  parallel  (mixture)  configurations  in  the 
electrode.  The  electrode  operation  is  considered  in  direct  current, 
and  the  capacitance  in  the  R//C  circuit  is  ignored  because  the 
relaxation  time  corresponding  to  the  R//C  circuit  is  negligible 
compared  to  the  duration  of  the  test  regimes  discussed  in  this  work. 
To  derive  resistance  values,  the  open  circuit  potential  (OCP)  versus 
state-of-charge  (SOC)  or  OCP  =  /( SOC)  curve  [1]  needs  to  be 
established.  This  curve  is  established  from  the  constant  current 
charge  and  discharge  regimes  at  a  very  low  rate,  such  as  C/25,  in 
a  half-cell.  The  electrode  resistance  value  versus  SOC  correspon¬ 
dence  is  derived  from  pseudo  Ohms’  law  approximation  at  any  rate 
C  and  SOC  as  follows: 


Rexp(C,SOC) 


l/ocp(SOC)-Uc(C,SOC) 

abs(/c) 


(1) 


where  Rexp(C,SOC)  is  the  electrode  resistance  as  a  function  of  rate  C 
and  SOC  as  derived  from  the  experimental  data.  Vocp(SOC)  and 
Vc(C,SOC)  are  electrode  potential  at  open  circuit  and  a  given  C  rate, 
respectively.  The  abs(Jc)  is  the  absolute  value  of  the  current  at 
a  given  C  rate.  Rexp(C,SOC)  is  the  sum  of  Ri  and  R2 ,  where  Ri 
corresponds  to  the  ohmic  resistance  calculated  from  the  initial 
voltage  drop  after  current  is  applied  [23]  and  R2  the  faradic  resis¬ 
tance  that  is  rate  and  SOC  dependent  in  the  ECMs  [22,23]. 

Flydro-Quebec  and  Timcal  kindly  provided  the  half-cell  data  for 
the  LFP  and  graphite  electrodes  at  different  rates,  respectively.  For 
the  LFP  electrode,  the  parameters  used  to  construct  the  sub-layer 
module  are  compiled  from  the  discharge  data  at  C/8,  C/4,  C/2,  C/1, 
2C  and  4C.  For  the  graphite  electrode,  the  parameters  are  compiled 
from  the  charge  data  at  C/50,  C/5  and  5C.  The  electrode  behavior  at 
any  other  rates  is  simulated  with  interpolations  or  extrapolations. 

Fig.  2  presents  the  simulation  of  the  electrode  behavior  for  (a) 
LFP  discharge  curves  and  (b)  graphite  charge  curves  at  rates  from  C/ 
25  to  2C  with  0.1  C  intervals.  For  the  convenience  of  discussion,  the 
potential  in  the  charge  discharge  curves  simulated  by  the  ECM  will 
be  denoted  as  Veiectrode@C/n,  where  C/n  is  the  imposed  rate  on  the 
half  cell.  The  model  can  simulate  electrode  behavior  with  high 
fidelity  to  extract  capacity,  resistance  and  kinetic  information 
accurately  for  any  rate  within  the  range  [23].  The  fidelity  can  be 
assessed  by  a  simple  verification,  as  shown  in  Fig.  3.  In  this  process, 
a  specific  rate  (e.g.  C/2  and  2C,  respectively)  is  removed  from  the 
parameter  matrix  for  the  sub-layer  module  construction.  The 
module  thus  constructed  shall  simulate  the  discharge  curve  for  this 
rate  and  the  result  is  compared  with  the  experimental  data.  As 
shown  in  the  figure,  respectively,  the  simulation  and  experimental 
data  for  C/2  and  2C  are  in  excellent  agreement. 

The  separation  of  current  density  at  the  electrode  and  cell  level 
is  of  great  importance  for  this  study  in  two  accounts:  (1 )  the  current 
density  applied  to  the  cell  should  be  treated  differently  in  the 
electrodes,  depending  on  the  apparent  area  used  by  each  electrode. 
(2)  The  current  density  at  each  electrode  may  also  vary,  depending 
on  the  aging  condition  and  degradation  mode.  Although  similar 
considerations  might  have  been  considered  in  the  electrochemical 
models  previously  published  by  others  [13-18],  it  has  never  been 
clear  to  us  how  such  an  aspect  was  incorporated  in  the  models.  We 
need  to  point  out  that  in  approaches  with  purely  experimental 
fitting,  it  is  quite  resource  demanding  to  consider  this  aspect 


a 


> 

<D 

O) 

£ 

o 

> 


4 

3.5 

3 

2.5 

2 


0  20  40  60  80  100 


SOC(%) 


Fig.  2.  (a)  Simulation  of  the  LFP  electrode  discharge  regimes  in  a  variety  of  rates  from 
C/25  to  4C.  (b)  Simulation  of  the  graphite  electrode  charge  regimes  from  C/25  to  4C. 


because  half-cell  data  on  a  single  electrode  need  to  be  tested  over 
an  extended  range  of  conditions  with  sufficient  details. 


2.2.  Full  cell  model 

To  emulate  the  cell  behavior  from  the  electrode  sub-layers 
effectively,  two  parameters  are  essential:  the  loading  ratio  (LR) 
between  the  negative  and  the  positive  electrode  (NE  and  PE);  i.e. 
LR  =  QneIQpe,  and  the  initial  irreversible  capacity  loss  of  the  NE  that 
compensates  the  SEI  formation. 

Fig.  4  presents  an  example  using  an  hypothetical  LFP  cell,  where 
an  original  LR,  LRq  =  1.45,  is  chosen  to  emulate  a  typical  high  power 
(TIP)  cell  design  [27].  Because  of  the  uncertainty  in  the  initial 


Fig.  3.  Simulation  of  C/2  and  2C  discharge  curves  of  a  Li  ||  LFP  cell.  The  data  of  each 
simulated  rate  (C/2  or  2C)  was  not  used  in  the  parameter  matrix  for  the  model 
construction  and  simulation,  but  as  the  basis  for  comparison  with  simulated  results. 
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Fig.  4.  Simulation  of  cell  voltage  as  a  function  of  SOC  with  the  formation  of  the  initial 
SEI  layer. 


Fig.  5.  Simulation  results  for  cell  discharge  regimes  at  C/25,  C/5,  C/2,  C,  2C  and  4C  using 
a  graphite  ||  LFP  cell. 


capacity  loss  at  the  NE,  it  is  more  convenient  and  reliable  to  use  the 
SOC  scale  at  the  PE  (SOCpe)  as  the  basis  for  the  calculation  of  cell 
SOC.  Thus,  the  Vne  =/( SOCne)  curve  can  be  transformed  and  scaled 
as  a  function  of  SOCpe,  Vne  =  /(SOCpe),  as  follows: 

SOCNE  *  (100%  -  SOCNE)LR0  (2) 

In  this  representation,  100%  SOCne  (fully  delithiated  NE)  corre¬ 
sponds  to  0%  SOCpe,  whereas  0%  SOCne  (fully  lithiated  NE)  should 
correspond  to  145%  SOCpe  because  of  the  LRo  =  1.45  in  the  capacity 
conversion.  All  SOC  values  are  presented  in  percentage. 

It  is  well  known  that  during  the  first  cycle,  the  NE  usually  loses 
some  capacity  due  to  the  formation  of  the  SEI  layer  at  the  surface 
[28].  In  Fig.  4,  an  irreversible  8.5%  active  material  loss  on  the  NE  is 
assumed,  which  is  equivalent  to  12%  capacity  loss  in  the  full  cell, 
based  on  the  LR  calculation  (i.e.  8.5%  x  1.45  =  12%)  and  the  fact  that 
the  full  cell  capacity  is  provided  by  the  lithium  content  in  the  PE. 
For  convenience  in  the  simulation,  an  SOCpe  offset  (OFS)  and  an 
apparent  LR  are  used  in  the  model  to  calculate  the  full  cell  behavior 
upon  initial  SEI  layer  formation.  This  initial  offset  (OFSini) 
corresponds  to  the  shift  of  the  VNe  in  correspondence  to  the  VPE  on 
the  SOCpe  scale  and  the  amount  of  graphite  loss  in  the  first 
cycle.  The  LR  accounts  for  the  graphite  loss  in  the  adjustment  of  the 
LRo,  i.e.: 


100%  —  OFSin 

Too% 


LRo 


(3) 


In  Fig.  4,  an  offset  of  OFSini  =  12%  SOC  and  LRjni  =  1.3  shall  yield 
approximately  the  same  composition  after  the  initial  irreversible 
capacity  loss  during  cell  formation  (from  8.5%  loss  in  NE  and 
LRo  =  1.45). 

Because  of  the  SOC  offset,  Equation  (2)  needs  to  be  modified  as: 


SOCne  =  (100%  -  SOCNE)LRini  +  OFSini  (2') 

Equation  (2')  is  generally  true,  including  situations  with  degra¬ 
dation  as  well. 

Fig.  5  presents  the  simulated  discharge  curves  of  the  full  cell 
described  in  Fig.  2  at  C/25,  C/5,  C/2, 1C,  2C,  and  4C.  The  discharge 
curves,  along  with  the  rate  capability  (as  expressed  by  Peukert 
curve)  and  the  polarization  resistance  (both  derived  but  not 
shown),  are  similar  to  those  measured  experimentally  for  cells  in 
a  commercial  HP  design  tested  in  our  laboratory  [27].  The  polari¬ 
zation  resistance  comprises  both  the  ohmic  resistance,  including 
contact  resistance  and  resistance  derived  from  the  conductance  in 
the  electron  conductive  network  of  the  electrodes  and  the  ionic 
conduction  in  the  electrolyte;  and  the  faradic  resistance  attributed 


from  charge  transfer,  mass  transport  in  the  porous  electrode 
matrices  and  ion  diffusion  in  the  active  materials. 

The  architecture  and  high  fidelity  of  the  model  allows  us  to  infer 
end-of-discharge/charge  (EOD/EOC)  SOCs  and  rest  cell  voltages 
(RCVs)  for  any  given  cut-off  conditions.  This  is  achieved  by  (1) 
calculating  the  pseudo-OCV  of  the  cell  (ps-OCVceii)  from  the  ps- 
OCPpe  or  ne  =  /( SOCpe)  curve  for  each  electrode,  which  is  obtained 
from  averaging  the  charge  and  discharge  curves  of  each  electrode 
measured  at  a  very  slow  rate  (such  as  C/25  or  lower)  [6];  (2) 
properly  matching  the  SOC  offset  and  LR;  thus,  the  full  cell  ps- 
OCVceii  =  /(SOCpe)  curve  is  constructed;  and,  (3)  converting  the 
SOCpe  to  an  SOC  scale  for  the  cell  according  to  Equation  (4): 


SOC^) 


SOCpE(Vceu)- SOCpe  (ocvggp) 

SOCpe  (ocvjgf)  -  SOCpe  (ocv^p) 10CU 


(4) 


The  EOD  SOC  depends  on  the  given  cutoff  condition  used;  e.g. 
2.5  V  in  this  study.  The  RCVeod  is  inferred  from  the  full  cell  ps- 
OCVceii  at  the  EOD  SOC. 


3.  Simulation  of  cell  degradation  modes 

Many  possible  cell  degradation  modes  have  been  reported  in  the 
literature  for  various  cell  chemistries.  It  appears  quite  common  in 
the  literature  that  the  loss  of  lithium  inventory  (LLI)  due  to  the 
formation  of  SEI  layers  is  often  observed  in  the  early  stage  of  cell 
aging  [4,6,14,16,17,20,29,30].  In  addition,  LLI  can  be  induced  by  Li 
plating  [31].  Another  common  degradation  mode  is  the  loss  of 
active  material  (LAM)  [4,6,14,16,17,20,29,30],  which  can  be  further 
categorized  into  four  types,  depending  on  the  affected  electrode 
and  the  degree  of  lithiation  (i.e.  predominantly  in  a  lithiated  or  de¬ 
lithiated  state)  in  which  the  LAM  occurs.  The  origin  of  LAM  could  be 
due  to  isolation  of  the  grains  of  the  active  materials  from  ionic  or 
electronic  conduction  network  in  the  electrode  [32],  dissolution  of 
transition  metal  into  the  electrolyte  solution,  change  in  the  elec¬ 
trode  composition  [33,34],  or  changes  in  the  crystal  structure  of  the 
active  material  (e.g.  due  to  decomposition  and  oxygen  evolution) 
that  result  in  a  reduction  in  the  extent  of  lithiation  [35].  Additional 
modes  for  capacity  loss  include  the  retardation  in  the  charge- 
transfer  kinetics  at  the  electrode  surface  [6,30,34,36],  ionic 
conduction  in  the  active  materials,  or  mass  transport  in  the  porous 
electrode  matrices,  which  attributes  to  the  increase  in  the  faradic 
resistance;  and  increase  of  the  ohmic  resistance  due  to  degradation 
in  electrode  contact  or  electrolyte  conduction,  as  commonly 
observed  in  polarization  resistance  increase  [13],  and  the  formation 
of  parasitic  phases  [35]. 
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By  properly  adjusting  parameters  in  both  layers  of  the  model,  all 
these  degradation  modes  can  be  simulated  with  sufficient  flexi¬ 
bility  to  accommodate  their  variations.  The  full  cell  layer  can  handle 
the  four  types  of  LAM  and  LLI,  whereas  the  half-cell  sub-layer  can 
handle  ohmic  resistance  increase,  faradic  rate  degradation,  and 
formation  of  parasitic  phases  that  are  electroactive. 


3.1  Loss  of  active  material  (LAM)  and  loss  of  lithium  inventory  (LLI) 

Fig.  6  presents  the  simulation  of  the  four  types  of  LAM  at  C/25  up 
to  50%  loss  of  capacity  (with  10%  increments,  in  dotted  lines).  The 
initial  and  the  state  of  50%  capacity  loss,  respectively,  are  denoted 
by  solid  lines.  Each  type  of  LAM  is  noted  by  the  electrode  in  which  it 
occurs  (in  the  subscript),  with  a  prefix  that  specifies  the  state  of 
lithiation;  i.e.  “li”  for  lithiated  and  “de”  for  delithiated. 

A  loss  of  active  material  (AM)  will  affect  the  LR.  Disregarding  if 
the  AM  is  lithiated  or  not,  as  long  as  some  PE  AM  is  lost,  the  LR  will 
increase.  Similarly,  the  LR  will  decrease  if  NE  AM  is  lost.  Given  a  50% 
AM  loss,  the  LR  will  double  if  some  PE  AM  is  lost  and  reduce  by  half 
if  some  NE  AM  is  lost.  The  evolution  of  the  LR  upon  a  given 
percentage  of  LAM  can  be  expressed  as: 


/100%  -  %LAMdeNE  -  %LAMliNE\ 
V 100%  -  %LAMdePE  -  %LAMliPE  j 


The  LAM  will  also  change  the  rate  on  the  affected  electrode, 
since  it  decreases  both  the  material  content  that  can  accept  charge 
and  the  active  surface  area  and  thus  increases  the  current  density. 
Therefore,  the  electrode-specific  rates  CpE  and  Cne  can  be  expressed 
as  a  function  of  the  “imposed”  rate  (Qmp)  on  the  cell,  the  extent  of 
degradation  (in  percentage)  in  the  electrode  and,  on  the  NE,  the  LR: 
100%Cimp 

CpE  =  100%  -  %LAMdePE  -  %LAMiiPE  (6) 


_  ioo%cimp 

NE  -  LR(100%  -  %LAMdeNE  -  %LAMliNE) 


(7) 


The  evolution  of  LR  and  the  CpE  and  Cne  are  not  affected  by  the 
state  of  lithiation  in  the  LAM.  The  electrochemical  behavior  of  the 
cell  Vceii  =/( SOCPE)  is  however  affected. 


In  the  charge  regime,  were  the  PE  suffering  from  LAM  (LAMhpe), 
less  charge  is  released  from  the  PE  when  it  reaches  EOC.  Conse¬ 
quently,  the  NE  should  have  a  higher  EOC  SOCne  than  that  of  the 
prior  cycle.  Vice  versa,  were  LAM  occurring  in  the  NE  (LAMdeNE)»  by 
the  time  the  PE  is  fully  delithiated,  the  NE  will  have  a  lower  SOCne 
than  that  of  the  prior  cycle,  since  a  less  amount  of  NE  will  accept  the 
same  charge.  Schematics  of  both  scenarios  are  presented  in  Fig.  6(a) 
for  the  LAMhpe  and  Fig.  6(c)  for  the  LAMdeNE,  respectively.  It  has  to 
be  noted  that  when  0%  SOCne  is  reached  earlier  than  100%  SOCpe  on 
the  PE  scale  (dotted  line  on  Fig.  6(c))  the  cell  is  at  risk  of  Li  plating 
(cf.  Section  3.2.3)  because  the  NE  is  fully  lithiated  before  the  PE 
completely  releases  all  its  Li  ions.  In  the  simulation,  as  the  condition 
of  EOD  does  not  change,  the  OFSini  remains  the  same  and  a  simple 
adjustment  of  the  LR  is  sufficient  to  emulate  both  cases. 

In  the  discharge  regime,  the  LAM-affected  electrode  will 
complete  its  electrochemical  reactions  at  a  higher  SOC.  Contrary  to 
the  charge  regime;  the  LAM  does  not  affect  the  state  of  the  affected 
electrode  in  the  beginning  of  the  discharge  regime  (i.e.  EOC  SOC). 
Fig.  6(b)  and  (d)  shows  the  case  for  the  LAMdep£  and  LAMhne, 
respectively.  Additional  adjustment  to  the  model  is  needed  to 
simulate  these  two  cases,  since  the  EOD  of  the  affected  electrode  is 
evolving  with  aging;  so  is  the  LR  matching  with  the  opposite 
electrode.  The  simulation  requires  adjustments  in  both  LR  and  SOC 
offset.  By  keeping  the  EOC  on  the  SOCpe  scale  unchanged,  we 
simulate  these  LAMs  by  adjusting  the  LR  and  SOC  offset  according 
to  the  percentage  of  LAM  and  the  relative  loading  matching  of  the 
two  electrodes  as  follows: 


1  R 

OFS  =  OFSini  +  LR  x  %LAMliNE  +  x  %LAMdePE  (8) 

bikini 

The  model  can  also  simulate  the  scenarios  involving  LLI.  LLI 
occurs  when  reversible  Li  ions  are  partially  consumed  by  parasitic 
reactions  in  the  cell.  LLI  commonly  occurs  during  charging  when  Li 
ions  are  used  to  grow  the  passivation  layer  on  the  NE.  It  may  also 
occur  during  discharging  if  a  passivation  layer  consuming  Li  is 
growing  on  the  PE.  Fig.  7  illustrates  the  case  for  a  charge  regime, 
where  LLI  does  not  affect  the  EOC  SOCpe,  since  the  PE  is  releasing 
the  Li-ions.  The  SOCpe  at  EOC  of  the  NE  is  however  evolving  with 
aging  and  becoming  lower  cycle  by  cycle;  since  less  Li  ions  are 
inserting  in  the  NE  each  time.  This  results  in  the  evolution  of  the 


Fig.  6.  Simulation  of  the  four  types  of  loss  of  active  material  (LAMs)  at  C/25  in  delithiated  PE  (dePE),  lithiated  PE  (liPE),  delithiated  NE  (deNE)  and  lithiated  NE  (liNE). 


M.  Dubarry  et  al.  /  Journal  of  Power  Sources  219  (2012)  204-216 


209 


Fig.  7.  Simulation  of  the  LLI  degradation  mode  at  C/25. 


loading  correspondence  in  the  discharge  curves  between  the  two 
electrodes,  as  shown  in  Fig.  7.  The  correspondence  can  be  adjusted 
by  including  the  percentage  of  LLI  in  the  calculation  of  the  SOC 
offset  as  follows: 

1 R 

OFS  =  OFSini +LR  x  %LAMliNE  +— —  x  %LAMdePE +%LLIch  -  %LLIdis 


The  subscripts  of  “ch”  and  “dis”  are  for  “charge”  and  “discharge” 
regimes,  respectively.  In  the  follow-on  discussion,  only  LLICh  is 
considered;  thus,  the  notation  LLI  refers  to  LLICh,  unless  specified. 


3.2.  Ohmic  resistance  increase ,  faradic  rate  degradation ,  and 
formation  of  parasitic  phase 

3.2.1.  Ohmic  resistance  increase  (ORI) 

Fig.  8(a)  shows  a  simulation  of  a  scenario  of  ohmic  resistance 
increase  (ORI)  by  100%  in  the  PE  and  66%  in  the  NE  and  their 
impacts  on  a  2C  discharge  curve  of  a  cell.  In  this  case,  the 
potential  of  each  electrode  is  altered  by  the  amount  of  ORI  with 
no  change  in  the  LR  or  the  SOC  offset.  Since  the  half-cell  ECM  in 
the  sub-layer  computes  the  initial  electrode  potential  response  at 
a  given  rate  C,  a  potential  correction,  due  to  the  additional  IR  drop 
induced  by  the  percentage  of  ORI  (%ORI),  needs  to  be  introduced 
as  follows: 

l/Cell@Qmp  =  VpE@CpE  -  CPERPE  ^  -  VnE@QmE 

+  Cne«ne(^)  (9) 


Fig.  8.  Simulation  of  (a)  ohmic  resistance  increase  (ORI)  at  2C  and  (b)  faradic  rate 
degradation  (FRD). 

becomes  more  gradual  than  before.  Such  consequences  make  the 
shape  of  the  discharge  curve  resemble  to  that  of  a  higher  rate  (at  2C 
in  the  example),  although  the  initial  ohmic  IR  drop  remains  the 
same.  Thus,  the  FRD  is  simulated  by  “downgrading”  the  reaction 
kinetics  to  a  lower  rate  (Crd)  than  the  initial  one  imposed  on  the  PE 
(CPE),  while  maintaining  the  initial  IR  drop  to  its  CPE  value,  and  by 
introducing  a  rate  degradation  factor  (RDF)  defined  as: 

RDFpE  =  ^2  (io) 

CpE 

in  the  rate  calculation  for  each  electrode  as  follows: 

£  _  100%RDFPECjmp  . 

PE  -  100%  -  %LAMdePE  -  %LAMliPE  1  j 


where  RPE  is  the  ohmic  resistance  associated  with  the  PE  and  RNE 
the  one  associated  with  the  NE  in  the  corresponding  Ri  in  the  sub¬ 
layer  module  ECMs,  respectively. 

3.2.2.  Faradic  rate  degradation  (FRD) 

Faradic  rate  degradation  (FRD)  occurs  when  an  electrode  does 
not  react  with  Li  ions  at  the  same  kinetic  rate  in  aging.  FRD  induces 
changes  in  the  electrochemical  behavior  of  the  cell,  such  as  IC  peak 
broadening,  and  increase  in  R2  in  the  corresponding  sub-layer 
module  ECMs.  Fig.  8(b)  shows  a  simulation  for  FRD  in  a  C/5 
discharge  of  a  LFP  electrode  due  to  aging.  Because  of  FRD,  the 
voltage  plateau  is  shortened  (i.e.  as  a  result  of  the  slower  kinetics,  so 
the  disparity  in  the  extent  of  reaction  throughout  the  electrode 
increases  with  aging)  and  the  voltage  decrease  toward  the  EOD 


^  _  100%RDFNEQmp  , 

NE  -  LRini(100%  -  %LAMdeNE  -  %LAMliNE)  1  } 

To  hold  the  initial  IR  drop  unchanged,  Equation  (9)  is  modified 
as  follows: 

Vcell @Cimp  =  VPE@CPE  -  CpeRpe  +  rdfpe  -  l) 

-Une@Cne  +  CneRne^  ioo%E  +  rc*^NE  ~  (^) 

The  same  approach  can  be  applied  to  the  temperature  effect  on 
reaction  kinetics.  As  the  temperature  of  an  electrode  changes,  the 
reaction  kinetics  often  vary  as  well. 


210 


M.  Dubarry  et  al.  /  Journal  of  Power  Sources  219  (2012)  204-216 


3.2.3.  Formation  of  parasitic  phase  (FPP) 

A  parasitic  phase  can  be  either  Li  reactive  or  not.  Either  case  can 
be  simulated  using  this  model.  If  a  non-reactive  parasitic  phase  is 
formed,  the  degradation  mode  can  be  modeled  as  a  LAM  on  the 
affected  electrode.  The  type  of  LAM  (on  lithiated  or  delithiated  AM) 
is  determined  by  the  mass  and  charge  balance  between  the  AM  and 
the  parasitic  phase  composition. 

If  the  parasitic  phase  is  Li  reactive,  the  affected  electrode  will 
resemble  a  “composite”  electrode,  showing  the  potential  changes 
of  all  electrochemically  active  phases  in  the  electrode  reactions.  The 
potential  and  the  associated  capacity  variation  in  the  electrode 
reaction  are  referred  to  as  “signature”  to  facilitate  following 
discussions.  The  signature  of  the  main  component  in  the  electrode 
should  diminish  gradually  through  the  course  of  degradation,  while 
the  signature  of  the  parasitic  phase  continues  to  grow  in 
a  progressively  apparent  fashion.  This  evolution  of  signature  can  be 
emulated  by  adding  sub-layers  in  the  model;  each  representing  the 
behavior  of  a  component  in  the  electrode  in  series  or  parallel  with 
one  another.  The  potential  of  the  “composite”  electrode  is  thus 
simulated  in  the  degradation  mode. 

An  example  is  provided  in  Fig.  9  for  the  LFP  electrode.  It  has  been 
reported  in  the  literature  that  exposing  LFP  raw  materials  in  the  air 
could  result  in  the  formation  of  LixFeP04(0H)x  phases  [37,38].  These 
phases  and  their  electrochemical  behavior  have  been  characterized 
in  Refs.  [39-41].  The  literature  data  of  a  LFP  cell  with  AM  aged  in 
the  air  can  be  adopted  in  the  model  and  the  cell  performance 
simulated.  The  data  for  the  hydrated  LFP  phase  in  Fig.  9  are  taken 
from  Ref.  [39]  where  the  C/50  charge  discharge  data  are  available. 
Fig.  9  presents  the  evolution  of  cell  voltage  vs.  SOCpe  curves  in 
a  graphite  ||  LFP  cell  with  0%  (thick  solid  line)  10%,  20%  (dashed 
lines)  and  30%  (thin  solid  line)  hydrated  phase  presence.  The 
hydrated  phase  is  electroactive  at  low  potentials;  thus,  the  shape  of 
the  PE  potential  curve  at  low  SOCs  is  modified  to  reflect  the  nature 
of  the  electrochemical  behavior  of  the  PE  composition  changes.  The 
hydrated  phase  can  only  take  up  to  0.8  Li  per  formula  even  at  low 
rates  [39].  Thus,  if  30%  of  the  hydrated  phase  is  present,  the  EOD 
SOC  should  be  offset  by  (l-0.8)*30%  =  6%.  The  presence  of  10% 
hydrated  phase  is  not  going  to  create  any  noticeable  changes  in 
capacity  with  the  existing  cutoff,  because  the  induced  SOC  offset 
and  changes  in  the  OCV  =  /(SOC)  curve  occur  in  the  range  of  OFSini. 
However,  if  30%  of  the  hydrated  phase  is  present,  it  will  induce 
a  noticeable  cell  signature  change  and  a  small  capacity  loss. 

Similar  phenomena  could  also  happen  in  the  NE;  for  instance,  Li 
plating  is  a  well-known  phenomenon.  Li  plating  often  occurs  under 
two  conditions:  (1)  if  the  IRne  drop  of  the  NE  during  charging 
makes  the  NE  potential  goes  <  0  V  or  (2)  if  the  LAMdeNE  is  signifi¬ 
cant  to  reach  a  condition  where  overcharging  in  NE  happens 


Fig.  9.  Simulation  of  a  graphite  ||  LFP  cell  discharge  curves  with  corresponding  PE 
potential  versus  SOCPE  curves  in  the  presence  of  0%  (thick  solid  line),  10%,  20%  (both  in 
dashed  lines),  and  30%  (thin  solid  line)  of  hydrated  phase,  respectively. 


(dotted  line  in  Fig.  6(c)).  If  the  plated  Li  is  passivated  progressively 
[42],  it  will  result  in  additional  LLI.  Therefore,  increasing  the  %LLICh 
in  Equation  (8')  shall  enable  the  simulation  of  this  scenario. 


3.3.  Multiple  processes  in  a  degradation  mode 

In  the  previous  sections,  the  simulation  of  individual  degrada¬ 
tion  mode  (e.g.  LAM,  LLI,  ORI,  FRD  and  FPP)  was  discussed.  By 
combining  all  the  equations,  this  model  is  capable  of  simulating 
degradation  modes  with  multiple  concurrent  processes  in  an 
electrode  or  both.  Fig.  10  presents  some  case  studies.  Fig.  10(a) 
shows  a  scenario  of  degradation  mode  comprising  a  logarithmic 
increase  of  LAMhne*  an  exponential  increase  of  LAMdePE  and  a  linear 
progression  of  LLI  over  100  cycles  of  aging.  The  model  computes  the 
evolution  of  the  contributions  to  capacity  fade  from  these  three 
phenomena  as  a  function  of  cycle  number  (Fig.  10(a))  and  simulates 
the  electrode  potential  and  cell  voltage  vs.  SOC  changes  from  the 
initial  to  the  condition  at  the  end  of  100-cycle  aging  (Fig.  10(b))  with 
proper  LR  and  SOC  offset  adjustments. 

This  versatility  to  simulate  complicated  degradation  modes  at 
any  rate  and  the  capability  of  predicting  capacity  loss,  EOD  SOC, 


Fig.  10.  (a)  Simulation  of  capacity  loss  versus  cycle  number  in  multiple  concurrent 
degradation  modes  of  LAMdepE,  LAMiiNE,  and  LLI,  and  (b)  the  corresponding  electrode 
potential  at  PE  and  NE  and  cell  voltage  versus  SOCPE  at  the  initial  and  final  stage  of 
degradation. 
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RCVeod,  (dQ/dV)  and  DV  (dV/dQ)  signatures  are  presented  below 
to  show  the  applicability  of  this  mechanistic  modeling  approach  for 
battery  diagnosis  and  prognosis. 

4.  Applications  to  battery  diagnosis  and  prognosis 

4.1.  Capacity  and  EOD  SOC  changes 

To  exemplify  the  utility  of  this  model  for  evaluation  of 
commercial  cells,  two  distinct  cell  designs  of  the  same  chemistry 
based  on  graphite  NE  and  LFP  PE  were  compared:  a  high  power 
(HP)  design  with  a  high  graphite  content  for  high  rate  capability 
and  a  high  energy  (HE)  design  of  a  lower  LR  with  low  rate  capability 
and  high  energy  content  (see  details  in  Ref.  [27]).  The  interest  in 
this  study  is  to  study  the  effect  of  LLI  on  capacity  retention  and  the 
corresponding  change  in  the  SOC  range  used  in  the  course  of  aging. 
This  analysis  is  motivated  by  our  previous  study  [4]  of  an  HE  cell  in 
which  the  cell  exhibited  a  unique  behavior  in  its  first  100  cycles  of 
aging  where  capacity  fading  was  observed  at  lower  rates  but  not 
higher  ones,  yet  the  cell  seemed  discharged  to  a  lower  SOC 
progressively.  It  is  curious  for  us  to  figure  how  the  LLI  results  in  such 
a  phenomenon  and  if  we  can  apply  the  knowledge  learned  from 
this  exercise  to  HP  design. 

Fig.  11(a)  presents  the  capacity  retentions  and  (b)  EOD  SOC 
evolution  simulated  for  50%  LLI  in  three  scenarios:  the  C/25 
discharge  curve  for  this  chemistry  and  the  2C  discharge  curves  for 
the  HE  and  HP  design,  respectively.  The  simulation  suggests  that 
the  trend  of  capacity  fade  for  the  HP  cell  @2C  is  similar  to  that  of  C/ 
25;  and,  it  is  proportional  to  the  amount  of  LLI.  In  the  HE  design,  the 
trend  of  capacity  fade  is  different  from  the  other  two.  At  first,  the 
capacity  @2C  is  much  lower  than  the  other  two  (because  of  the  HE 
design)  but  stays  stable  over  a  range  of  cycles,  although  LLI  is 
recurring  in  the  cell.  Upon  reaching  15%  LLI,  the  capacity  fade 
begins  to  appear  and  follows  the  same  trend  as  the  other  two.  The 
simulation  results  are  consistent  with  our  experimental  observa¬ 
tions  reported  in  Ref.  [4].  Fig.  11(b)  shows  the  evolution  of  the  EOD 
SOC,  which  tends  to  increase  slightly  for  C/25  and  HP@2C 


discharges  but  decreased  by  20%  for  the  HE@2C  discharge  before 
stabilizing  and  following  the  same  trend  with  the  other  two. 

The  peculiar  behavior  of  the  HE  design  can  be  explained  by 
plotting  the  evolution  of  the  electrodes  V = /(SOC)  in  all  three  cases. 
Fig.  11(c)  and  (d)  present  the  evolution  of  the  potential  of  PE,  NE  and 
the  cell  voltage  in  the  course  of  degradation  for  the  HP  and  HE 
designs,  respectively.  The  dotted  lines  represent  the  initial 
OCV  =  /(SOC)  curves  for  both  electrodes.  The  solid  lines  are  the 
initial  discharge  curve  of  the  cell  at  2C  and  the  corresponding 
potential  curves  for  the  electrodes.  The  dashed  lines  are  the 
discharge  curves  showing  the  capacity  fade  induced  by  consecutive 
LLI  in  10%  increments  up  to  50%  of  LLI. 

In  the  HP  design  (Fig.  11(c)),  the  cell  can  deliver  capacity  close  to 
its  maximum  at  2C  and  the  capacity-limiting  electrode  is  the  NE. 
This  is  not  the  case  for  the  HE  design.  Because  of  the  lower  inhered 
rate  capability,  the  capacity  is  initially  limited  by  the  PE. 

As  described  in  Section  3.1,  due  to  LLI  the  NE  loading  progres¬ 
sively  shifts  toward  higher  SOCs.  Since  the  NE  is  the  limiting  elec¬ 
trode  for  the  HP  design,  its  capacity  fade  will  be  reflected 
immediately  in  the  cell  capacity  loss;  on  the  contrary  in  the  HE 
design,  as  long  as  the  PE  is  limiting  the  capacity,  the  NE  capacity 
fade  would  not  be  detected.  The  capacity  fade  will  appear  only 
when  the  NE  loading  is  shifted  sufficiently  to  make  the  NE  capacity 
limiting. 

The  evolution  of  the  EOD  SOCs  in  the  aging  process  can  be 
explained  similarly.  For  the  HP  design,  the  cutoff  condition  at  EOD  is 
determined  at  the  NE,  and  this  scenario  remains  the  same  in  the 
aging.  Thus,  the  EOD  SOC  for  the  C/25  discharge  is  almost  the  same, 
although  gradually  sliding  to  higher  values  is  predicted.  The  sliding 
of  SOC  toward  high  values  for  both  C/25  and  2C  discharges  is 
imputable  to  the  overall  capacity  loss  and  an  unaltered  rate  capa¬ 
bility.  The  unaltered  rate  capability  suggests  the  same  difference  in 
capacity  at  the  EOD  between  C/25  and  2C  through  aging.  As  the 
maximum  capacity  is  decreasing  in  the  course  of  aging,  this 
constant  capacity  difference  represents  a  progressively  higher  EOD 
SOC.  For  the  HE  design,  the  capacity  difference  between  C/25  and 
2C  is  not  constant  but  decreasing  through  the  course  of  aging. 


Fig.  11.  (a)  Capacity  retention  and  (b)  EOD  SOC  evolution  as  a  function  of  %LLI  and  electrode  potentials,  cell  voltage  as  a  function  of  SOC  simulated  with  a  hypothetical  50%  LLI  for  (c) 
HP  and  (d)  HE  designs. 
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Therefore,  the  EOD  SOC  at  2C  decreases  until  the  NE  becomes  the 
capacity-limiting  electrode. 

4.2.  Cell  diagnosis  and  prognosis 

Since  the  model  can  produce  results  for  inference  of  battery 
degradation  modes,  it  is  a  powerful  technique  for  battery  diagnosis 
and  prognosis.  Here  the  derivation  of  two  useful  pieces  of  infor¬ 
mation,  the  incremental  capacity  (IC)  and  differential  voltage  (DV) 
curves,  is  illustrated  to  show  the  utility  of  this  approach.  Both  IC  and 
DV  signatures  have  been  used  as  inference  techniques  capable  of 
revealing  battery  degradation  mechanisms  and  quantifying  the 
degree  of  degradation  in  a  cell.  The  ability  to  derive  these  signa¬ 
tures  is  an  advancement  in  the  inference  techniques,  since  degra¬ 
dation  results  can  now  be  simulated  under  various  scenarios  to 
help  us  identify  the  mechanisms  with  much  less  uncertainty.  For 
instance,  in  the  comparison  of  the  four  types  of  LAMs,  we  can 
identify  which  mode  might  occur  in  the  cell  using  dQJdV  or  dV/dQ. 
analysis,  without  any  ambiguity.  In  some  cases  where  some 
degradation  modes  may  give  similar  results  in  the  capacity  loss,  we 
can  simulate  each  hypothesis  in  the  model  to  see  if  subtle  differ¬ 
ences  can  be  identified  among  various  scenarios.  Thus,  in  principle, 
through  careful  analysis,  additional  testing,  and  simulation,  we 
should  be  able  to  test  various  hypotheses  and  identify  the  mecha¬ 
nism  with  more  quantitative  information;  and,  verify  and  enhance 
the  inference  technique. 

A  cell  voltage  comprises  a  proper  matching  of  the  Vne  =  /(SOC) 
and  VpE  =  /(SOC)  curves  based  on  their  loading  ratio  correspon¬ 
dence  and  degradation  modes,  as  executed  in  the  top  layer  module 
of  the  model.  The  evolution  of  the  cell  voltage  in  an  aging  process  is 
thus  attainable  from  the  simulation.  Fig.  12  presents  the  C/25 
discharge  curves  with  30%  capacity  loss  for  the  four  LAM  and  LLI 
modes.  These  degradation  modes  not  only  show  a  variety  of 
capacity  loss  scenarios,  from  almost  none  in  30%  LAMdeNE  to  about 
50%  in  30%  LAMune,  but  also  produce  a  variety  of  VCeii  =  /(SOC) 
signatures.  These  variations  in  cell  performance  characteristics  are 
the  basis  for  battery  diagnosis  and  prognosis  by  inference. 

Fig.  13  shows  the  variations  in  the  IC  signature  for  six  degra¬ 
dation  modes:  i.e.  30%  in  four  LAMs  and  LLI,  and  100%  ORI.  In  the 
figure,  the  heavy  solid  line  represents  the  initial  state,  the  thin  solid 
line  the  final  state,  and  the  dotted  lines  1  /5th  increments  in  the 
degree  of  degradation.  Fig.  14  shows  the  same  results  for  the  DV 
analysis.  To  facilitate  discussions,  the  same  convention  as  in  our 
previous  work  will  be  used  [5,6].  The  graphite  staging  phenomena 
are  numbered  as  0  to  0,  where  0  is  the  reaction  from  LiCi2  to  LiC6 
and  so  on,  and  0  the  last  staging  transition  leading  to  the  deli- 
thiated  graphite  (as  shown  in  Fig.  12).  In  the  dQJdV  =  /(V)  curve, 
these  five  staging  phenomena  coupled  with  the  LFP  potential 
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Fig.  12.  V  =  f(Q)  for  an  HP  cell  at  the  initial  state  and  with  30%  of  degradation  from 
various  modes  comprising  the  four  types  of  LAMs  and  LLI  at  C/25. 


plateau  are  observed  as  five  IC  peaks  [4].  In  the  DV  analysis,  these 
five  staging  phenomena  become  valleys  in  the  dV/dQ=f(Q )  curve. 
There  are  four  dV/d Q=f(Q)  peaks  representing  non-stoichiometry 
in  the  single-phase  regions  (solid  solution),  as  noted  by  A  to  D.  The 
peak  correspondence  is  labeled  in  the  upper  left  plot  (a)  in  each 
figure. 

The  utility  of  these  IC  and  DV  curves  is  further  explained  here. 
As  shown  in  the  left  column  of  Figs.  13  and  14,  LAMhpe,  LAMdeNE 
and  LLI  show  similar  IC  and  DV  signatures.  They  share  a  common 
theme  that  the  graphite  cannot  be  lithiated  to  the  same  level  as  it 
was  initially;  so,  the  main  feature  in  the  respective  signatures  is  the 
loss  of  intensity  in  the  first  peak  0  on  the  IC  curve  and  the  valley 
0  on  the  DV  curve.  For  the  LAMhpe  and  LLI,  the  intensity  of  peak/ 
valley  0  is  affected  only  when  0  is  completely  lost.  In  the  LAM¬ 
deNE,  all  ©  to  0  peaks/valleys  are  fading  from  the  beginning,  thus 
LAMdeNE  can  be  unambiguously  identified.  It  is  interesting  to  note 
that  it  is  almost  impossible  to  distinguish  between  LAMhpe  and  LLI 
using  the  DV  analysis  in  this  case.  It  is  however  possible  in  the  IC 
analysis,  since  peak  0  is  shifting  toward  higher  voltages  faster  for 
LLI  than  for  LAMhpe,  although  it  might  be  difficult  to  decipher 
between  the  two  from  the  C/25  IC  signature,  especially  for  small 
capacity  fade  (<5%)  and  if  the  PE  potential  is  really  a  constant 
potential  plateau. 

There  is  a  shift  in  the  trend  of  signature  evolution  in  LAMdePE 
when  the  PE  becomes  the  capacity-limiting  electrode.  Before  the 
shift,  no  changed  are  observed;  beyond,  the  last  graphite  staging 
process  cannot  be  completed  anymore;  thus,  the  peak  0  in  the  IC 
curve  and  peak  D  in  the  DV  curve  are  disappearing.  This  disap¬ 
pearance  is  followed  by  higher  losses  inpeaks/valleys  ©  and  Q. 

In  the  case  of  LAMdeNE,  all  the  IC  peaks  are  shrunk  (since  less 
graphite  is  involved  in  the  reactions)  except  peak  0  that  increases. 
This  is  because  LR  >  1,  and  with  LAMdeNE,  more  and  more  capacity 
in  the  staging  0  plateau  is  accessible  in  the  intercalation.  In  the  DV 
representation,  the  accessible  capacity  in  valley  0  increases;  thus, 
all  the  other  valleys  are  shrinking,  because  the  total  capacity  is 
constant.  Hence,  all  DV  peaks  are  moving  toward  lower  capacities. 

The  lower  right-hand  plots  in  Figs.  13  and  14  are  the  results  of 
ORI,  which  indicate  a  shift  of  all  IC  peaks  toward  lower  voltages  on 
the  IC  curves;  in  contrast,  no  change  at  all  for  DV  peaks  in  the  DV 
curves.  This  is  due  to  the  fact  that  the  ORI  only  affects  the  cell 
voltage,  not  the  capacity.  Therefore,  without  capacity  fade,  the 
effect  can  only  be  seen  on  a  voltage  scale. 

In  summary,  both  IC  and  DV  analyses  offer  merits  and  draw¬ 
backs  for  battery  degradation  inference.  One  should  consider  their 
merits  to  choose  proper  tools  to  highlight  different  aspects  of  the 
degradation  phenomena.  Although  capacity  quantification  might 
be  easier  with  the  DV  analysis,  IC  analysis  may  inhere  a  better 
sensitivity  to  decipher  degradation  modes.  It  reflects  cell  degra¬ 
dation  signature  on  a  voltage  scale,  which  provides  a  better  refer¬ 
ence  to  the  state  of  the  battery  than  the  capacity  scale  that  varies 
with  aging. 

This  study  shows  that,  in  the  case  of  graphite  1 1  LFP  chemistry, 
the  derivation  of  signature  for  various  cell  degradation  modes  is 
rather  feasible,  but  it  is  difficult  to  distinguish  between  LAMhpe  and 
LLI  unambiguously.  This  is  due  to  the  fact  that  the  LFP  used  in  the 
example  has  a  flat  potential  plateau  and  good  rate  capability.  We 
trust  that  the  two  modes  would  be  distinguishable  in  other 
chemistries,  or  if  the  LFP  electrode  is  carefully  characterized  with 
better  resolution  in  rate  and  temperature  effects. 

5.  Implications  to  battery  testing 

In  addition  to  the  derivation  of  IC  and  DV  signatures  and  their 
evolutions  in  aging,  the  model  can  allow  us  emulate  other  inter¬ 
esting  features  to  aid  the  understanding  of  battery  testing  data 
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Fig.  13.  IC  signatures  from  the  initial  state  (thick  solid  line)  and  at  the  end  of  30%  degradation  (thin  solid  line)  with  1  /5th  increment  (dotted  lines)  for  (a-d)  each  of  the  four  LAMs, 
(e)  LLI,  and  (f)  a  100%  ohmic  resistance  increase  (ORI)  at  C/25. 


usually  received  from  typical  cell  evaluation  procedures  [1—6]. 
Following  are  some  useful  illustrations  to  explain  such  aspects: 


5.1.  OCV  =  f  (SOC)  evolution 

Fig.  15(a)  presents  the  evolution  of  ps-OCV  =  /( SOC)  curves 
simulated  for  different  LAMs  and  LLI,  where  Fig.  15(b)  shows  the 
evolution  of  the  ps-OCV  =  /(SOC)  curve  with  30%  LLI  in  3%  incre¬ 
ments.  The  examples  reveal  that  there  are  variations  in  the  ps- 
OCV  =  /(SOC)  curves,  depending  on  the  modes  and  the  degree  of 
degradation.  As  being  revealed  in  the  test  data  [6],  degradation 
modes  may  change  the  correspondence  between  the  PE  and  NE 
loading. 

Understanding  such  cell  voltage  evolution  in  aging  is  important 
because  one  of  the  most  effective  way  to  determine  the  cell  SOC  is 
to  measure  the  rest  cell  voltage  (RCV)  and  use  an  OCV  =  /(SOC) 
curve  to  infer  the  SOC  value  [6,43].  If  the  correlation  between  the 
RCV  and  SOC  varies  with  aging,  the  inference  of  SOC  should 
gradually  lose  its  accuracy.  This  issue  thus  introduces  a  very 
difficult  hurdle  to  overcome  in  battery  control  and  management. 
For  instance,  in  Fig.  15(a),  it  is  shown  that  for  a  RCV  =  3  V,  the  SOC 
of  the  cell  is  7%  at  the  initial  state  but  could  be  in  the  range  of  2% 
and  7%  after  aging  depending  on  degradation  scenarios.  This 
indicates  that  to  accurately  determine  the  SOC  of  the  cell  from  the 
RCV,  it  is  essential  to  calibrate  the  ps-OCV  =  /(SOC)  as  often  as 
possible  through  reference  performance  tests  (RPT)  at  low  rates 
upon  aging  [6],  which  is  quite  impractical.  Another  example 
showing  the  issue  with  the  evolving  EOD  SOC  =  /(RCVEod)  is 
presented  in  Fig.  15(c),  in  which  the  RCV  evolution  is  shown  with 


up  to  30%  LLI  in  a  FIE  cell  (following  the  case  of  Fig.  11(b)). 
Although  the  EOD  SOC  evolution  seems  straightforward  with 
a  steady  decrease  over  aging,  the  RCVeod  evolution  is  affected  not 
only  by  the  voltage  plateaus  (see  Q  in  Fig.  15(c))  where  the  OCV  is 
constant  over  a  range  of  SOCs  but  also  by  the  evolution  of  the  ps- 
OCV  =  /(SOC)  curves.  Indeed,  the  RCV  of  2.9  V,  which  should 
correspond  to  2.5%  SOC  initially,  now  is  5%  SOC  after  aging.  These 
evolutions  of  both  RCVs  and  SOCs  have  been  observed  experi¬ 
mentally  for  commercial  LFP  cells  [4]. 

Our  model  could  be  used  as  a  calibration  tool  to  calculate  the  ps- 
OCV  =  /(SOC)  curves  evolution  at  any  stage  of  degradation,  which 
can  be  validated  between  two  RPTs,  for  diagnosis.  It  could  also  be 
used  for  prognostics  by  extrapolation  from  previous  RPTs,  if  the 
degradation  modes  were  identified  from  the  test  data  and  the 
results  quantified  to  feed  the  model.  Since  our  model  does  not 
require  much  computation,  the  results  can  be  easily  used  to 
compose  lookup  tables  (from  single-electrode  half-cell  data);  thus, 
our  approach  can  be  easily  implemented  into  battery  management 
systems  (BMS)  to  improve  SOC  and  SOFI  tracking. 


5.2.  Effect  of  the  electrode-specific  rate  on  the  cell  behavior 

Another  interesting  capability  derived  from  this  study  is  the 
ability  to  assess  the  effect  of  “actual”  rate  on  an  electrode  with 
LAMs,  especially  for  cells  that  have  low  rate  capability  (as  in  the  FIE 
designs).  It  should  be  noted  that,  in  conventional  test  regimes,  the 
rate  is  calculated  from  the  rated  capacity  of  the  cell  and  imposed  at 
the  cell  level.  It  is  often  not  attended  at  the  electrode  level  (as  the 
“actual”  rate  on  the  electrodes).  This  disparity  in  the  rate  could  lead 
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Fig.  14.  DV  signatures  from  the  initial  state  (thick  solid  line)  and  at  the  end  of  30%  degradation  (thin  solid  line)  with  1  /5th  increment  (dotted  lines)  for  (a— d)  each  of  the  four  LAMs, 
(e)  LLI,  and  (f)  a  100%  ohmic  resistance  increase  (ORI)  at  C/25. 


to  significant  impacts  on  the  electrode  performance,  while  it  might 
be  less  noticeable  at  the  cell  level.  Fig.  16  presents  an  example 
where  the  evolution  of  the  IC  signatures  with  30%  LAMupe  are 
simulated.  In  the  figure,  the  heavy  solid  line  represents  the  initial 
state,  the  thin  one  the  final  state,  and  the  dotted  ones  results  in  6% 
increments.  It  is  clear  that  the  IC  peaks  are  slowly  shifting  toward 
lower  voltages.  This  would  be  interpreted  as  a  polarization  resis¬ 
tance  increase,  if  a  real  cell  were  tested,  whereas  no  ORI  or  FRD 
were  included  in  the  simulation.  The  shift  of  the  IC  peaks  is  in  fact 
related  to  the  increase  of  the  “actual”  rate  on  the  PE  that  goes  from 
about  2C  to  3C  (considering  the  same  “imposed”  rate  on  the  cell, 
whereas  only  2/3  of  the  active  material  on  the  PE  are  reacting  to  the 
current).  Indeed,  as  shown  in  Fig.  5,  when  the  “imposed”  rate  on  the 
cell  is  increased,  the  IR  drop  shall  cause  the  IC  peaks  shift  toward 
lower  voltages.  This  observation  implies  that  the  polarization 
resistance  increase  often  reported  in  the  literature  upon  degrada¬ 
tion  may  not  be  entirely  due  to  SEI  layer  growth  or  degradation  in 
charge-transfer  kinetics;  to  the  contrary,  some  may  be  imputable  to 
LAM  that  changed  the  “actual”  rate  at  the  electrode  that  is  induced 
by  degradation. 

Interestingly,  the  same  phenomenon  could  be  interpreted  for 
kinetic  degradation  as  often  observed  in  aging  [6,30,36],  since  the 
EOD  SOC  is  also  increased  if  the  rate  is  increased  (Fig.  5).  More¬ 
over,  LLI  could  also  be  responsible  for  an  increase  of  EOD  SOC 
toward  the  end  of  life  (Fig.  11(b)).  In  this  case,  the  RCVEod  should 
be  stable,  since  the  SOC  evolution  is  caused  by  the  reduction  of 
the  maximum  capacity  in  the  cell,  not  a  change  in  kinetic 
performance. 


Further  investigations  are  in  progress  to  de-convolute  the  effect 
of  LAM  on  the  FRD  and  ORI  using  some  test  data  from  commercial 
cell  aging  experiments. 

5.3.  Second  stage  of  cell  degradation 

In  the  testing  of  commercial  cells,  it  is  often  noticed  that  the 
capacity  fade  is  rather  constant  in  the  first  stage  with  a  degradation 
mode  associated  with  LLI.  An  accelerated  capacity  fade  is  observed 
in  the  second  stage  of  degradation;  typically,  above  500  cycles  [4,6] 
in  today’s  cell  designs.  The  second  stage  is  often  assumed  to  come 
from  kinetic  degradation  accompanied  with  LAM  [4,6,17,30,36].  The 
hypothesis  mentioned  in  Section  5.2  could  explain  the  co¬ 
appearance  of  the  LAM  and  kinetic  changes,  but  it  could  not 
explain  why  they  appear  so  suddenly. 

Figs.  6  and  7  provide  some  explanations  for  the  sudden 
appearance  of  the  second  stage  degradation.  Some  of  the  LAMs 
may  be  masked  for  up  to  30%  of  capacity  fade,  if  they  occur  on  the 
electrode  that  is  not  capacity  limiting  to  the  cell  performance. 
However,  such  degradation  will  eventually  catch  up  with  that  of 
the  limiting  electrode  and  shift  the  roles  in  the  capacity  limiting 
mechanism  over  aging;  thus,  these  LAMs  finally  become  notice¬ 
able.  Worse  than  that,  since  the  LAM  might  be  significant  in 
magnitude  already,  its  impact  on  the  cell  performance  could 
appear  more  severely,  as  in  an  accelerated  manner.  Our  model 
simulation  suggests  that  a  degradation  mode  may  start  benign  but 
progress  in  a  more  aggressive  manner  in  the  aging;  thus,  it 
becomes  prominent  after  certain  incubation.  This  might  be  the 
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Fig.  15.  Evolution  of  the  ps-OCV  =/( SOC)  curve  for  (a)  30%  capacity  fade  from  each  of 
the  four  types  of  LAM  and  LLI,  (b)  30%  HI  with  3%  increments,  and  (c)  associated 
evolution  in  RCV  and  EOD  SOC  as  a  function  of  %LLI. 


case  for  the  LAMdePE  that  incubates  for  a  while,  when  LLI  is 
occurring  in  a  cell.  Since  LLI  shifts  the  NE-to-PE  loading  corre¬ 
spondence  in  the  same  direction  as  LAMdePE  (as  shown  in  Figs.  6 
and  7),  the  impact  on  capacity  fade  from  LAMdePE  might  be 


Fig.  16.  IC  signatures  simulated  with  30%  LAMiiPE  for  the  HE  design  with  2C  aging.  The 
heavy  solid  line  represents  the  initial  state,  the  thin  one  the  final  state,  and  the  dotted 
ones  results  in  6%  increments. 


Fig.  17.  Evolution  of  the  capacity  fade  in  a  cell  with  linear  LLI  and  exponential  loss  of 
LAMdePE- 

masked  until  it  becomes  prominent.  Our  experience  suggests  that 
LLI  impacts  capacity  fade  rather  linearly  in  aging,  whereas  LAM 
may  follow  a  power-law  or  exponential  dependence  [8].  This 
scenario  is  simulated  in  Fig.  17  where  a  linear  0.03%/cycle  LLI  is 
accompanied  by  an  exponential  LAMdePE-  After  400  cycles,  the 
capacity  fade  is  still  driven  primarily  by  LLI  despite  >  25%  LAMdePE- 
After  500  cycles,  the  LAMdePE  starts  to  dominate  the  capacity 
fade  at  a  higher  pace.  This  capacity  fade  scenario  is  similar  to 
the  aging  behavior  of  commercial  cells  observed  in  the  laboratory 
tests  [4,6]. 

6.  Conclusion 

For  practical  applications,  one  should  realize  that  only  a  few 
parameters  (likely  time-resolved  voltage,  current,  and  tempera¬ 
ture)  might  be  available  from  the  cell  monitoring  in  the  battery 
management  systems  (BMS),  considering  all  costs,  logistics,  and 
scheduling  issues  in  on-board,  in-line  battery  diagnosis  and  prog¬ 
nosis.  Therefore,  any  invasive  characterizations  or  uses  of  compli¬ 
cated  instrumentation  are  unrealistic  and  cost-prohibitive.  It  is 
essential  to  design  simple  diagnostic/prognostic  tools  and 
approaches  for  on-board,  in-line  battery  monitoring  and  control 
that  do  not  require  costly,  complicated  monitoring/sensoring 
devices  and  computation  overhead.  To  understand  degradation 
modes  in  commercial  cells  over  aging,  a  unique  simulation  capa¬ 
bility  is  demonstrated  here.  This  modeling  approach  is  built  from 
individual  electrode  behavior,  including  degradation  modes,  and 
proper  adjustment  of  the  loading  ratio  between  the  two  electrodes 
in  a  cell  design  to  achieve  high  fidelity  in  the  simulation.  This 
simulation  capability  offers  great  values  to  allow  battery  diagnosis 
and  prognosis  in  BMS.  The  advantages  of  this  approach  are: 

■  A  unique  inference  technique  for  battery  diagnosis  and 
prognosis, 

■  A  mechanistic  model  that  can  synthesize  a  variety  of  cell  aging 
scenarios  based  on  degradation  modes,  including  loss  of  active 
material,  loss  of  lithium  inventory,  kinetic  degradation  or 
increase  of  polarization  resistance,  formation  of  parasitic 
phases,  Li  plating,  and  any  combination  of  them, 

■  A  physicochemical  process-based  model  that  constitutes  from 
electrode-specific  degradation  modes  to  construct  cell  perfor¬ 
mance  and  degradation  scenarios  with  consideration  of  cell 
designs  and  operating  conditions, 

■  A  versatile  modeling  tool  without  constraints  on  chemistry 
variations,  cell  designs,  battery  sizes  and  geometries,  and 
operating  or  aging  conditions. 
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Because  of  this  capability,  we  illustrated  some  interesting 
degradation  modes  that  were  difficult  to  distinguish  in  the  exper¬ 
iments  but  can  be  easily  deciphered  in  this  approach.  For  instance, 
loss  of  active  materials  could  occur  without  detectable  signatures  in 
the  beginning  of  aging  but  surface  after  an  incubation  period,  as 
observed  in  some  experiments  with  commercial  cells. 
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